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Abstract 

The spatially anisotropic triangular antiferromagnet is investigated with the numerical diagonal- 
ization method. As the anisotropy varies, the model changes into a variety of systems such as the 
one-dimensional, triangular, and square-lattice antiferromagnets. Taking into account such a geo- 
metrical character, we impose the screw-boundary condition, which interpolates smoothly the one- 
and two-dimensional lattice structures. Diagonalizing the finite clusters with N = 16, 20, . . . , 32 
spins, we observe an intermediate phase between the VBS and Neel phases. Suppressing the in- 
termediate phase by applying the ring exchange, we realize a direct VBS-Neel transition. The 
simulation data indicate that the transition is a continuous one with the correlation-length critical 
exponent v = 0.80(15). These features are in agreement with the deconfinement-criticality scenario 
advocated by Senthil and coworkers in the context of the high-temperature superconductivity. 
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I. INTRODUCTION 



According to the deconfinement-criticality scenario j-^^ in (2 + 1) dimensions, the phase 
transition separating the VBS and Neel phases is continuous, accompanied with uncon- 
ventional critical indices. Naively, such a transition should be discontinuous,- because the 
adjacent phases possess distinctive order parameters such as the VBS-coverage pattern, and 
the sublattice magnetization, respectively. A good deal of field-theoretical investigations 
have been made to clarify this issue.-^ For instance, as a low-energy effective description, 
the QED 3 theory has been investigated;-^ it would be intriguing that the theory exhibits a 
deconfinement transition.— On the one hand, because of the magnetic frustration, the Monte 
Carlo simulation suffers from the negative-sign problem to realize the VBS phase. However, 
in recent Monte Carlo simulations,—*^ the biquadratic interaction, rather than the magnetic 
frustration, has been utilized. Thereby, it turned out that the biquadratic-interaction-driven 
transition is a continuous one with unconventional critical indices. (On the contrary, in Refs. 
000 , it was claimed that the transition would be a weak first-order one.) 

In this paper, we investigate the spatially anisotropic triangular antiferromagnelr^ 1 ^ with 
the ring exchange by means of the numerical diagonalization method. As the anisotropy 
varies, the model changes into the one- and two-dimensional systems, and correspondingly, 
a variety of phases appear. To cope with such a geometrical peculiarity, we impose the 
screw-boundary condition, which interpolates the one- and two-dimensional lattice struc- 
tures smoothly. 

To be specific, we present the Hamiltonian for the spatially anisotropic triangular anti- 
ferromagnet with the ring exchange; 

h = j$> • s i + f s i ■ s ; + J 4 E( p ^< + p ^)- C 1 ) 

(ij) ((ij)) mi] 

The quantum spin-1/2 operators {Si} are placed at each triangular-lattice point i. The sym- 
bol Pijki denotes a ring-exchange operator with respect to a plaquette k, 1} consisting of 
two adjacent triangles; namely, as to a plaquette state the operation P1234 trans- 

lates it into Isl^sl)- The summations Yl(ij)> S«y))> anc ^ Stofc/l run over an possible vertical 
nearest-neighbor pairs, remaining nearest-neighbor pairs, and plaquette spins, respectively; 
the triangular lattice is directed so that one of the triangular edges points upward. The 
parameters J, J', and J4 are the corresponding coupling constants. (In the next section, we 
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present an explicit expression for the Hamiltonian matrix, referring to the technical details 
of the screw-boundary condition.) Hereafter, we consider J as a unit of energy; namely, we 
set J = 1. 

In Fig. [Tj we present a schematic phase diagram; the details are explained in Sec. IIHI 
As mentioned above, the aim of this paper is to survey the direct VBS-Neel transition; 
in this sense, the ring exchange J4 is significant to realize the VBS-Neel transition. A 
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23l : First, the case J4 



number of limiting cases were studied in Refs. 
was investigated with the spin-wave,— series-expansion,— >22 large-A,— and numerical- 
diagonalization^ methods. The regime of the intermediate (triangular antiferromagnetic) 
phase 24 was estimated as 0.27 < J' < 2, 0.25 < J' < 1.43, 0.13 < J' < 1.71, and 0.78(5) < 
J' < 1.15(10), respectively. (Some analyses predict two types of intermediate phases. Such 
a detail is ignored for simplicity.) These results appear to be unsettled. It is a purpose of 
this paper to survey the intermediate phase. Second, the spatially isotropic (J' = J) case 
in the presence of the ring exchange was investigated in Ref. [23J; here, the generic types of 
ring-exchange interactions were considered in the context of the Helium adsorbate. It was 
reported that the J4-driven phase transition occurs in agreement with our observation. 

As mentioned above, the model ([I]) has a geometrical peculiarity. That is, as the spatial 
anisotropy J' changes, the model ([1]) reduces to the one-dimensional (J 7 = 0), triangular 
(J 7 = 1), and square-lattice (J 7 — > 00) antiferromagnets successively. (Hence, for sufficiently 
large J' , the conventional non-collinear Neel phase appears.) Notably enough, the phase 
diagram, Fig. [U reflects this geometrical character. In order to take into account this 
geometrical character, we implemented the screw-boundary condition, which interpolates 
the one- and two-dimensional-lattice structures smoothly. 

In fairness, it has to be mentioned that the VBS-Neel transition was studied for the 
frustrated square-lattice antiferromagnet, namely, the J1-J2 model.— According to the 
series-expansion method,— the Neel {J%j J\ < 0.4), VBS (0.4 < J2/J1 < 0.6), and collinear 
(0.6 < J2/J1) phases appear successively, as the magnetic frustration changes. The VBS 
phase seems to be dominated by the presence of the collinear phase. (Note that for J2/J1 — > 
00, the system reduces to two independent square-lattice antiferromagnets. The collinear 
state consists of two independent Neel orders.) In this paper, we dwell on the triangular 
antiferromagnet flTJ), which exhibits an isolated VBS-Neel transition. 

The rest of this paper is organized as follows. In Sec. HIl we explicate the simulation 
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algorithm, placing an emphasis on the screw-boundary condition. In Sec. 1111} we show the 
finite-size-scaling analysis of the simulation data. In Sec. HVt we present the summary and 
discussions. 



II. SIMULATION METHOD: SCREW-BOUNDARY CONDITION 

In this section, we present an explicit expression for the Hamiltonian, Eq. (TO , under the 
screw-boundary condition. 

To begin with, we present a schematic drawing of the finite-size cluster in Fig. [2J As shown 
in the figure, the spins constitute a one- dimensional (d = 1) alignment {Si} (i = 1, 2, . . . , N). 
The dimensionality is lifted to d = 2 by the bridges over the vth-neighbor interactions. As 
mentioned in the Introduction, the spatially anisotropic triangular antiferromagnet possesses 
a geometrical character such that it reduces to a one-dimensional antiferromagnet in the limit 
J' — > 0. In this sense, the geometrical peculiarity is seized by the screw-boundary condition. 
Actually, for a rectangular cluster with the system size 6x6, for instance, the length of 
the independent chains in the limit J' — > is merely L = 6. On the contrary, owing to 
the screw-boundary condition, we attain treating the chain length L = 32 along the J-bond 
direction. 

To be specific, we present an explicit expression for the Hamiltonian matrix. We propose 
the following expression; 

H = JH(1) + J'(H(v) + H(v + 1)) + J 4 (# 4 (l, v) + #4(1, v + l) + H A (v + 1, v)). (2) 
Here, the wth-neighbor Heisenberg interaction H(v ) is given by 

N 

H(v) = J2Si-Si+v (3) 

i=l 

(The periodic condition, namely, S^+i = Sj, is imposed.) Similarly, the ring exchange is 
introduced via 

N 

H 4 (j, v) = ^{P i)i+j) i +V)i+j+v + h.c). (4) 



i=l 

We set the screw pitch to 



, nUN) + 1 for N > 24 
v(N)={ 1 ; " (5) 

n(y/~N) otherwise 



with the round-off function n(x) = [x + 0.5] and Gauss' symbol [...]; i.e., n(2.4) = 2. The 
screw pitch v(N) converges to v(N)/y/~N — > 1 for large system sizes N — > oo; hence, the 
spins form a VN x \fN network embedded on the torus. The rule, Eq. (EJ), is intended 
to suppress the finite-size errors; actually, by Eq. ([HJ), we can set the screw-pitch v to an 
even number (for small N), which turns out to improve the finite-size behavior even for 
small system sizes. More specifically, the screw-boundary condition introduces a frustration 
particularly for the Neel-type magnetism (J 7 — > oo), and the frustration effect is suppressed 
by the above rule, Eq. (JSJ). 

The above formulae complete the basis of our scheme. As shown in Fig. [2J the embedding 
geometry under the screw-boundary condition is essentially one-dimensional, admitting us to 
calculate the Hamiltonian- matrix elements systematically with Eq. (j2J). In the next section, 
utilizing the Lanczos method, we diagonalize the Hamiltonian matrix for the system sizes 
N < 32. 

III. NUMERICAL RESULTS 

In this section, we present the numerical results. We calculate the excitation gap 

AEi{k, S* tot ) = Ek(k, S z tot ) - E (0, + ), (6) 

with the ith low-lying energy E^k, S* ot ) (i — 1, 2, ... ) within the sector (k, S* ot ). Here, the 
index k denotes the wave number within the Brillouin zone — ir < k < tt. We impose the 
screw-boundary condition (Fig. [2]), and the Bloch wave k extends along the spiral (J-bond) 
chain; hence, the reciprocal space is one-dimensional. The quantum number Sf ot denotes 
an eigenvalue of the operator Yl!i=\ ^t- m ^ ne case °f ^tot = 0> additionally, we introduce 
an index ±, which specifies the inversion symmetry with respect to Sf — > —Sf. The sector 
(0, + ) contains the ground state. In this sector, we shift the % index so as to express the 
ground-state energy as E (0,0 + ) via i — > i — 1. (The ground-state energy is the starting 
point of all excitations, and it is sensible to index the ground-state energy as i£j=o rather 
than i = 1.) The linear dimension L of the cluster is given by 

L = VN, (7) 

because the N spins constitute a two-dimensional network as shown in Fig. [2J 
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A. Spatially anisotropic triangular antiferromagnet: J4 = 



In this section, we survey the regime without the ring ex ch ange J 4 = 0; as mentioned in 
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and the details of the 



the Introduction, this case has been studied in Refs. 
intermediate phase remain unclear. 

In Fig. [31 we plot the excitation gap AEi(n,0 + ) for J 4 = 0, various J' and N = 
16, 20, . . . , 32. We notice that the level crossings take place at J' « 0.65 and J' « 1.1. That 
is, the softening instability, AEi(n, + ) < 0, occurs in the intermediate regime. We estimate 
the range of the intermediate phase as 

0.65(15) < J' < 1.1(1). (8) 

Here, as an error indicator, we utilize the data scatters of the J'-intercept among N = 20, 
24, 28, and 32. (Several related studies are overviewed afterward.) 

Surveying various parameter ranges, we found that the elementary-excitation gap opens 
at either k = or n. The softening of the branch k = it suggests that the magnetic order 
along the J-bond direction is unstable against a staggered modulation. Such a staggered 
modulation fits the boundary condition (constraint) such that the chain length N is always 
set to an even number. On the one hand, as shown in Fig. [21 the number of spiral turns, N/v, 
of the chain is a fractional number, and the magnetism along the spiral direction may not fit 
the embedding geometry. Hence, the staggered order along the chain direction becomes even 
stabilized, resulting in the k = it softening. On the one hand, in the VBS phase J' < 0.65, 
the energy gap AEi(tt, + ) gets closed as the system size enlarges; eventually, the ground 
state may be doubly degenerated in the thermodynamic limit. This double degeneracy 
suggests that the J-bond chain is covered by the dimers. (In this sense, the VBS picture of 
the present system is not so complicated, as compared to that of the square lattice.— 1 ^) 

On the contrary, in the Neel phase 1.1 < J', a positive gap AEx(tt, + ) > starts to 
open. In fact, in the limit J' —>■ 00, the model reduces to the square- lattice antiferromagnet. 
Hence, the spins along the diagonal (J-bond) direction align ferromagnetically, and the 
k = it excitation exhibits a mass gap. 

It is a good position to make an overview of the related studies. According to the 
spin-wave,— series-expansion,— large-A,— and diagonalization studies,— the range of the 
intermediate phase is estimated as 0.27 < J' < 2, 0.25 < J' < 1.43, 0.13 < J' < 1.71, 
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and 0.78(5) < J' < 1.15(10), respectively. Our result, Eq. (jSJ), indicates that the VBS 
phase persists up to a considerably large J', suggesting that the VBS phase is robust. 
Similar conclusion was drawn from the diagonalization study by Weng and coworkers.- 3 
They diagonalized the rectangular clusters with the sizes 6 x 4, 8 x 4, and 6x6. Such 
a rectangular geometry is suitable for investigating the Neel-type magnetic structure. On 
the contrary, the screw-boundary condition meets the quasi-one-dimensional system (VBS 
phase). The agreement between these approaches would be encouraging. 

As a reference, in Fig. HI we present the ground-state energy per unit cell, E o (0, + )/N, 
with N = 32 for the same parameter range as that of Fig. [3j In the small- J' regime, 
the ground-state energy is close to the Bethe-ansatz solution, E /N = —0.443 . . . , for the 
one-dimensional Heisenberg antiferro magnet. This fact suggests that the VBS phase is of 
one-dimensional character. 

Last, we mention a number of remarks concerning the phase diagram. We made similar 
analyses for various values of J4 ~ 0, The result is summarized in Fig. [TJ as suggested by Eq. 
(jHJ), the intermediate-phase boundaries are not determined very precisely, and the boundaries 
in Fig. [1] are only schematic. (The critical branch separating the VBS and Neel phases 
is considered in the next section.) Concerning the singularities of the phase boundaries 
surrounding the intermediate phase, we followed the claim 22 that the intermediate- VBS 
(-Neel) phase transition is discontinuous (continuous). 



B. Suppression of the intermediate phase by the ring exchange J4 

In this section, we explore the regime with the ring exchange J4 7^ 0. 

In Fig. [5], we plot the excitation gap AEi(n,0 + ) for J 4 = 0.07, various J' and N = 
16, 20, . . . , 32. Comparing the result with that of Fig. [21 we notice that the ring exchange J4 
suppresses the intermediate phase (softening instability). As mentioned in the Introduction, 
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at J' = 1. 



the suppression of the intermediate phase by J 4 was demonstrated in Ref. 
In the present study, we have yet another parameter J', and we are able to investigate the 
J'-driven VBS-Neel transition. 

In Fig. [6l we plot the scaled energy gap L z AEi(0,l) for J4 = 0.07, various J', and 
N = 16, 20, . . . , 32. Note that the sector S* ot = 1 corresponds to the triplet excitation 
created preferentially on the J-bond chain. The behavior of the triplet excitation contains 
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information on the phase transition from the VBS phase. Here, we set the dynamical critical 
exponent to z — 1, following the conclusion of the Monte Carlo analyses^i^ for the square- 
lattice antiferromagnet. (Note that the energy gap AE± is proportional to the reciprocal 
correlation length, and the product LAE± is a dimensionless (scale invariant) quantity.) 
According to the scaling theory, the curves of the scaled energy gap should intersect each 
other at the critical point. In fact, we observe that a phase transition takes place at J' = 1.2. 
Surveying various parameter ranges, we arrive at the phase diagram, as depicted in Fig. [TJ 
A remark is in order. As mentioned in the Introduction, naively, the VBS-Neel transition 
should be discontinuous;- actually, the adjacent phases possess distinctive order parameters 
such as the dimer-coverage pattern and the sublattice magnetization, respectively. However, 
according to the deconfmement-criticality scenario,- the spinon deconfming from a sea of 
singlet dimers changes the singularity to a continuous one. Our result supports this scenario. 
In the next section, surveying a critical index, we investigate the criticality of the VBS-Neel 
transition more in detail. 



C. Deconfinement criticality 

In this section, we estimate the critical exponent v for the VBS-Neel transition. 

In Fig. d we present the finite-size-scaling plot, (J' - J'JL^-LAE^Q, 1), for J 4 = 0.07, 
various J', and N = 16, 20, . . . , 32. Here, we set the scaling parameters to J' c = 1.2 and 
v = 0.8; note that the former parameter J' c = 1.2 was determined in Fig. [6j The data of Fig. 
[7] collapse into a scaling curve, confirming that the transition is indeed critical. Moreover, 
the critical exponent acquires an enhancement, as compared to that of the 3c? Heisenberg 
universality, v = 0.7112(5). 32 (An overview of the related studies is addressed afterward.) 
Similarly, in Fig. [HJ we present the finite-size- scaling plot, (J' - J'^L^-LAE^O, 1), for 
J' c = 1.4, v = 0.8, J4 = 0.1, and N = 16, 20, ... , 32. Again, the data collapse satisfactorily. 
Surveying various parameter ranges, we arrive at 

v = 0.80(15). (9) 

This is a good position to make an overview of the preceding Monte Carlo studies. As for 
the square-lattice antiferromagnet , the biquadratic-interaction-driven VBS-Neel transition 



was investigated in Refs. 
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and 
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and the critical exponent was estimated as v = 0.78(3) 



S 



and v = 0.68(4), respectively. Moreover, as for the quasi-one- dimensional spin-1 antiferro- 



31 



for a field-theoretical interpretation. 



magnet, the index v = 1/2.9 was reported;— see Ref. 
(Note that these models are free from the negative-sign problem, and the quantum Monte 
Carlo method is applicable.) We notice that the results are not quite settled. A notable point 
is that the exponent^ v = 0.78(3) is significantly larger than that of the d = 3 Heisenberg 
universality class v = 0.7112(5),— suggesting a peculiarity of the deconfinement criticality. 
Our result, Eq. (jUJ), also suggests a tendency of an enhancement as to v. Nevertheless, 
our simulation result provides an evidence that the VBS-Neel transition is a critical one in 
agreement with the deconfinement-criticality scenario advocated by Senthil and coworkers. 

IV. SUMMARY AND DISCUSSIONS 

The spatially anisotropic triangular antiferromagnet with the ring exchange, Eq. (pQ), was 
investigated by means of the numerical diagonalization method. As the spatial anisotropy J' 
varies, the model changes into a variety of systems such as the one-dimensional, triangular, 
and square-lattice antiferromagnets successively. Taking into account such a geometrical 
character, we adopt the screw-boundary condition, as shown in Fig. [2j 

First, we survey the regime without the ring exchange J4 = 0. The simulation result 
indicates that the intermediate phase appears in 0.65(15) < J' < 1.1(1). Our result shows 
that the VBS phase is robust^ 1 ^ against the interchain coupling J'. Second, by applying 
the ring exchange J 4 , we suppress the intermediate phase. Eventually, we attain the di- 
rect VBS-Neel transition, which is under the current theoretical interest in the context of 
the high-temperature superconductivity. Postulating z = 1,—^ we analyze the simulation 
data in terms of the finite-size scaling. Thereby, we estimate the correlation-length criti- 
cal exponent as v = 0.80(15), confirming that the VBS-Neel transition is indeed a critical 
one. The exponent is comparable to the preceding Monte Carlo results, v = 0.78(3)^ and 
v = 0.68(4),— calculated for the square-lattice antiferromagnet. 

Our result provides an evidence that the VBS-Neel transition is critical, realizing the 



23l . generic 



deconfinement criticality. Here, the ring exchange plays a significant role. In Ref. 
types of ring-exchange interactions are considered in the context of the Helium adsorbate. 
Such an extension may also lead to an improvement as to the finite-size behavior. This 
problem will be addressed in a future study. 
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22 M.Q. Weng, D.N. Sheng, Z.Y. Weng, and R.J. Bursill, Phys. Rev. B 74, 012407 (2006). 

23 G. Misguich, C. Lhuillier, B. Bernu, and C. Waldtmann, Phys. Rev. B 60, 1064 (1999). 

24 B. Bernu, P. Lecheminant, C. Lhuillier, and L. Pierre, Phys. Rev. B 50, 10048 (1994). 

25 J. Oitmaa and Zheng Weihong, Phys. Rev. B 54, 3022 (1996). 

26 J. Sirker, Z. Weihong, O.P Sushkov, and J. Oitmaa, Phys. Rev. B 73, 184420 (2006). 

27 D. Poilblanc, A. Lauchli, M. Mambrini, and F. Mila, Phys. Rev. B 73, 100403(R) (2006). 

28 P. Sindzingre, Phys. Rev B 69, 094418 (2004). 

29 M. Mambrini, A. Lauchli, D. Poilblanc, and F. Mila, Phys. Rev. B 74, 144422 (2006). 

30 K. Harada, N. Kawashima, and M. Troyer, J. Phys. Soc. Japan 76, 013703 (2007). 

31 F.S. Nogueira, S. Kragset, and A. Sudbo, Phys. Rev. B 76, 220403(R) (2007). 

32 M. Campostrini, M. Hasenbusch, A. Pelissetto, P. Rossi, and E. Vicari, Phys. Rev. B 65, 144520 
(2002). 



11 




FIG. 2: A schematic drawing of the spin cluster under the screw-boundary condition for the 
spatially anisotropic triangular antiferromagnet (pQ) is presented. As indicated, the spins constitute 
a one-dimensional (d = 1) alignment {Si} (i = 1,2, ... , N) via the longitudinal coupling J. The 
dimensionality is lifted to d = 2 by the bridges over the -uth-neighbor pairs through the transverse 
coupling J'. The screw pitch v is given by Eq. ([5]). 
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FIG. 3: The energy gap AE 1 (ir,0 + ), Eq. j6]), is plotted for J 4 = 0, various J', and N = (+) 16, (x) 
20, (*) 24, (□) 28, and (■) 32. A softening instability, AE\ < 0, takes place in the intermediate 
phase 0.65(15) < J' < 1.1(1). 
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FIG. 4: The ground-state energy per unit cell, Eq(0, + )/N, with N = 32 is presented for the same 
parameter range as that of Fig. [3l In the small- J' regime, the ground-state energy is close to the 
Bethe-ansatz solution, Eq/N = —0.443 . . . , for the one-dimensional Heisenberg antiferromagnet. 




FIG. 5: The energy gap AE^n, + ), Eq. ©, is plotted for J 4 = 0.07, various J', and N = (+) 16, 
(x) 20, (*) 24, (□) 28, and (■) 32. Owing to the ring exchange, the softening instability, namely, 
AE\ < 0, does not occur any more. 
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FIG. 6: The scaled triplet-excitation gap LA£i(0, 1) is plotted for J4 = 0.07, various J', and N = 
(+) 16, (x) 20, (*) 24, (□) 28, and (■) 32. A continuous transition takes place at J' c = 1.2. 




14 




15 



